SPECTRAL ANALYSIS OF NON-COMMUTATIVE HARMONIC 
OSCILLATORS: THE LOWEST EIGENVALUE AND NO CROSSING 

FUMIO HIROSHIMA AND ITARU SASAKI 

Abstract. The lowest eigenvalue of non-commutative harmonic oscillators Q{a, f3) 
r^' • (a > 0,/3 > 0, a/3 > 1) is studied. It is shown that Q(q:,/3) can be decomposed into 

four self-adjoint operators, 

5_l _ <T=±,p=l,2 

Qh' and all the eigenvalues of each operator Q^-p arc simple. We show that the lowest 

■^ I eigenvalue E of Q(a, /3) is simple whenever a ^ (3, and E is also the lowest eigenvalue 

1^^ ■ of Q+i (rasp. Q+2) when /? > a (resp. a > 13). Furthermore a Jacobi matrix 

^vq I representation of Qa-p is given and spectrum of Q^p is considered numerically. 
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1. Introduction 



J2 [ The non- commutative harmonic oscillator is introduced by A. Parmeggiani and 

M. Wakayama |PW01t IPW02t IPW03] as a non-commutative extension of harmonic 
oscillators. We also refer to the monumental book |ParlO] by A. Parmeggiani. It is 
defined by 

• ! as an operator in "H = C^ ® L^(]R). Here A, J & Mat2(K), J is positive definite 

Q I symmetric, and B skew-symmetric. Furthermore A+iJ is positive definite. It is shown 

^ ' in |PW02l IPW03J that A and J can be assumed to he A = (^ J^ J , J = ( J ~ M , 

and a and /3 satisfy 
|>^ ; a > 0, /3 > 0, a/3 > 1. (1.2) 

c^ ! We fix A and J as above, and throughout this paper we assume (II. 2p . Under f ll.2p . Q 

is self-adjoint on the domain D{Q) = C'^^{D{cP/dx^)r\D{x^)) and has purely discrete 
spectrum Eq < Ei < E2 < ■ ■ ■ /^ 00. When a = /3, Q{a, (3) is equivalent to the direct 
sum of a harmonic oscillator. Then Ej = -Ej+i for j = 0, 2, 4, ■ ■ ■ . On the other hand, 
when a ^ P, Q{a,P) is regarded as a g-deformation of harmonic oscillators, and the 
spectral analysis of Q is nontrivial. 

An eigenvector associated with the lowest eigenvalue E = Eq is called a ground 
state in this paper. A long-standing problem concerning eigenvalues of Q{a, (3) is to 
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determine their multiplicity explicitly. Let a^ (3. Let En = -E.„(a,/3) denote the n-th 
eigenvalue of Q{a,f3). The map c„ : {a, (3) i— )■ En{a,(3) G M is called an eigenvalue- 
curve. To consider the multiplicity of eigenvalues is reduced to studying crossing or 
no crossing of eigenvalue-curves. 

We state a short history concerning studies of the multiplicity of eigenvalues of Q. 
In |PW03] it is shown that the multiplicity of any eigenvalues of Q is at most three and 
an alternative proof is given in |Och01j . In a numerical level it is found in |NNW02] 
that eigenvalue-curves cross at some points but the lowest eigenvalue is simple. The 
multiplicity of eigenvalues of Q is also considered in |IW07j , where it is derived that 



1\ r .. /a/3-1 ^ r. f 1\ r ^. /a/3-1 

n- -] mm.{a,[3\^ — < E2n~i < E2n < "- - 7^ maxja,/?}^ 



for n = 1, 2, 3, ■ ■ ■ . From this we can see that the multiplicity of E is at most two if 
P < 3a OT a < 3/3. In |Par04] it is shown that E is simple but for sufficiently large 
o;/3. Furthermore in |HS12] it is proven that the lowest eigenvalue is at most two and 
all the ground state are even for (a,/3) G D^, where D^ = {(a,/3)|a,/3 > \/2}, and 
it is also shown that E is simple for {a, (3) & D with some subset D C D^. Recently 
Wakayama |Wakl2] breaks through in studying the multiplicity of E. It is proven 
that if all the ground states are even, then E is simple whenever a ^ (3. Combining 
|Wakl2j with |HS12] . it is immediate to see that E is simple for (a, (3) G D^. 

In this paper we settle down the question concerning the multiplicity of the lowest 
eigenvalue of Q, i.e., we prove that E is simple for all values of a and [3 [a ^ (3) in 
Theorem 13.11 Moreover no crossing between eigenvalue-curves associated with an odd 
eigenvector and an even eigenvector is also proven in Corollary 15.21 

This paper is organized as follows. In Section HI we decompose Q{a, (3) into four self- 
adjoint operators: Q{a, (3) = ^^^^ „^i 2 Qap- It is shown that each Q^-p is equivalent 

to some Jacobi matrix Qap, and all the eigenvalues of Q^p are simple. In Section |3l 
we show that the lowest eigenvalue of Q{a,/3) is simple. In Section HJ we construct 
a unitary transformation Q^p \-^ Q„p such that e~*'^'"p is positivity improving, and 
it is shown that the ground state is in a positive cone. In Section El we show that 
Q-p — Q+p > A(a,/3), p = 1,2, with some A{a,/3). In particular, if A{a,(3) > 0, 
then there is no crossing between the n-th eigenvalue- curve of Q-p and that of Q+p. 
In Section [6l we show some numerical results. 

2. Decomposition of Q{a,(3) and Jacobi matrix 

2.1. Decomposition of Q{a,/3). Let a = -js{x + ^) and a* = -^{x — ^) be the 
annihilation operator and the creation operators, respectively. In terms of a and a*, 
Q can be expressed as 

Q = A{a*a + -) + -{aa- a*a*). (2.1) 

Let T-L^ (resp. "H.) be the set of even (resp. odd) functions in "H, and P+ (resp. P_) 
be the orthogonal projection onto 71+ (resp. Ti-). Let |?7,) be the n-th normalized 
eigenvector of a*a, i.e., |n) = ^(a*)" |0) with |0) = 7r"^/^e~^ /^. Let C |n) be the 
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one-dimensional subspace spanned by \n) over C Hence the Wiener-Ito decomposition 
L'^iR) = ®'^^QC\n) follows. The total Hilbert space is 



n 



C |n)' 
C|n) 



n=0 J n=0 

We use this equivalence without noticing. Since a\n) = ^/n\n — l) and a* \n) = 
\/n + 1 |n + 1), we see that aa : Tin — >■ 'H„_2 and a*a* : Tin — ;■ 'Hn+2- Furthermore 
a* a leaves Tin invariant. Then we have Q : Tin -^ 'Hn-2 ® Hn ® 'Hn+2- From these 
observation we can find invariant domain of Q. We denote the orthogonal projection 
onto C In) by |n)(n|, and define orthogonal projections on "H by 

P,in) . (I"H"I «) . n(") ^ (» ,„»„,) . (2.2) 

Note that 1 = Xl^o(-^t(^) ©-^il'^))- ■'■^ order to decompose Q, we define the following 
orthogonal projections: 



^+1 = Er=o(^t(4r^) © Pii^n + 2)), T+2 = Er=o(n(4n) © P^{^n + 2)), 

^-1 = Er=o(^t(4r2 + 1) © P4(4r^ + 3)), T_2 = Er=o(n(4n + 1) © P^{^n + 3)). 



Since |2n) is even and |2n + 1) is odd, one has T+i + T+2 = P+ and T_i + T_2 = P_. 
We set "Ho-p = Ran(To-p). Then "H is decomposed as 

"H = (^ ^crp- (2-3) 

cr=±,p=l,2 

Theorem 2.1. Operator Q is reduced by Tiap, cr = ±, p = 1, 2. 

Proo/. We see that o^Pj{n) D Pj{n-2)a'^, a*a*Pj{n) D Pj{n + 2)a*a* and a*aPj(n) D 
Pj{n)a*a for all n = 0, 1, 2, ■ ■ ■ , and j =t, |. Clearly it holds that APj{n) = Pj{n)A, 
JP^{n) = Pi{n)J and JPi{n) = P^{n)J. Then QT^-p D T^-pQ and the theorem 

follows. n 

Let us set Qo-p = Q\Hap- Then it holds that 

Q= Q.p. (2.4) 

o-=±,p=l,2 

2.2. Jacobi matrix representation of Qap- We construct a unitary operator im- 
plementing equivalence between Q^^ and a Jacobi matrix. Set 

oo oo 

?7+i = X;(^t(M © niS'^ + 2)) - Y.^P^i^'' + 4) ® ^■^(S'^ + 6))- (2-5) 

n=0 n=0 

This operator is unitary on "H+i and we have 



Q+i = U~lQ+iU+i = T+i ('A(a*a + ^) - ^{aa + a^a*)") T+i, 



(2.6) 
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Figure 1. RanT+i is supported on "■" 
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Figure 2. RanT+2 is supported on "■" 


n 





1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 




t 


n 


■ 


n 


n 


n 


■ 


D 


n 


n 


■ 


n 


n 


n 




i 


n 


n 


n 


■ 


n 


n 


n 


■ 


n 


n 


n 


■ 


n 




Figure 3. RanT„i is supported on "■" 
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Figure 4. RanT_2 is supported on 



In a similar way to f/+i one can define the unitary operators 



t/+2, t/-i and t/_2 on 1-1+2, H-i and 71-^2, respectively, such that 

Q+2 = t/+2Q+if/+2 = T+2 ( A{a*a + ij) - i^i^'O. + a*a*) j T+2, 

g_i = UzlQ-iU-i = T_i (A{a*a + ^) - ^{aa + a*a*)\ T_i, 

Q_2 = U:^Q^2U-2 = r_2 (A{a*a + ^) - ^{aa + a*a*) j T_2. 

For sequences a = (ao, ai, a2, ■ ■ ■ ) and b = {bo,bi, ,b2,- ■ ■), we define the Jacobi matrix 

/bo ao \ 

flo bi ai 

ai 62 ■ • • 



J{a,b) 



\0 



(2.7) 



■••/ 



which acts in the set of square summable sequences, i'^. Set a^ = (ao-(O), acr(l), ■ ■ ■ ) 



and b, 



crp 



(6.p(0),fe.p(l), 



where 



a+{n) = -^/i2^{TT)i2^^^T2), a.{n) = - V(2n + 2)(2n + 3), n = 0, 1,2, ■■■ 

«(! + 4n) for even n 



6+1 (n) 
b-i{n) 



/3(1 + 4n) for odd n, 

a(3 + 4n) for even n 
/3(3 + An) for odd n. 



6+2 (n) = 6+1 (n) 
6_2(n) = 6_i(n) 



{a,/3)^(/3,«) 



(a,/3)^(/3,«) 
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For (T = ± and p = 1, 2, we define the Jacobi matrix Q^-p by 

1 



Q 



o-p 



-J{a„,k 



apj 



(2.8) 



Set 
Let & 



]4n) 




iKj'. 





|4n + 2) 

'T=o e ^■ 



>n 



0,1,2, 



is a complete orthonormal system of "H+i. 



be tlie standard basis of i'^. We define tlie unitary operator 



y+i : H+1 ^ £2 by r+i 



\An) 




e2„ and F+i 



e2n+i- Tlien one can 





\4n + 2) 

compute tlie matrix element of Q+i as Q+i = y+iQ+i^-^T/- Similarly one can define 
the unitary transformations such that the following theorem holds. 



Theorem 2.2 (Jacobi matrix representations). For a = ±,p = 1,2, the operators 
Qcrp are unitarily equivalent to the Jacobi matrix Qap- 

Remark. In the case oi a = /3, Q^i = Qa2 for a = ±. Explicitly each Qo-p is expressed 
as 



Q, 



/1-2 



Vl-2 

5/3 

V3^ 



0\ 



V5^ 13/3 

-VTs 



-V7-i 
17a 



/9-10 



lo 



/9-10 21^ 



Q 



_ 1 

-2-2 



/3 -v^r2 

Vr2 5a -V3^ 

-V3^ 9/3 -v^S^ 

-Vl^ 13a -VTS 

-VTS 17/3 -v^gTio 



0\ 







/9-10 21a 



Q-1 - 2 



3a -V2^ 

-v^2l? 7/3 -Vi^ 

-VF5 11a -V(r7 

-v^fTY 15/3 -v^8^ 



0\ 



/F9 19a -VlO-11 







VUMT 23^ 



Q- 



313 -V2li 
-\/2li la -yJTl) 

-VF5 11/3 -Ve^ 
-V6^ 15a 



0\ 



/8^ 19/3 -VICTTT 



-VlO-ll 23a 



VO . ./ 

Theorem 2.3. Each eigenvalue of Q„p, a = ±, p = 1,2, is simple. 
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u..)^_n e f. Then A 



Proof. Let A be any eigenvalue of Q+i with an eigenvector u — v"'n;„=o 
and u satisfy the recurrence relations: 

Un+i = a+{n)~^ {{X - b+i{n))un - a+{n - l)un^i) , neNo, (2.9) 

M_i = 0. (2.10) 

Note that a_|_(n) ^ 0. Solutions of system (I2.9p - fl2.10l) are uniquely determined by the 

term uq G C. Hence the multiplicity of any eigenvalue of Q+i is simple. Proofs for 

other cases are similar. D 



3. Simplicity of the lowest eigenvalue of Q{a,(3) 

In this section, we state the main theorem in this paper. 

Theorem 3.1. Assume that a ^ /3. Then the lowest eigenvalue of Q{a,l3) is simple 
and the ground state is even. 



In order to show Theorem 13.11 we introduce a remarkable proposition given by 
Wakayama |Wakl2] . 

Proposition 3.2. Assume that (1) a ^ /3; (2) all the ground states of Q{a,/3) are 
even, i.e., ker((5(a,/3) — E) G H^. Then the lowest eigenvalue of Q{a,P) is simple. 

Let Qa = Qai ® Qa2, fT = +, — . Then Q is decomposed into the direct sum of even 
part and odd part, Q = Q+ © Q-. Let E^j = mi Spec{Qa). 

Lemma 3.3. The inequality E^ < E^ holds. 



Proof Let $_ 



$-1 

$-2 



be a normalized ground state of Q-. Note that ^-a, a = 1, 2, 
are odd functions. We define even functions $_ G ?/+ by 

^-^{x), ifx>0. 

Note that $_p G D{{-d'^ /dx'^f/^) and 



$_ 



$-1 

$-2 



^-Ax) 



if X < 0. 



||(ci/rfx)<l_f 
Thus one has 



\\{d/dx)^^ 



'^^^''^'dx^^P 



<l'_p/,X — $. 

dx 



P,P 



1,2. 



(3.1) 



E+ < (g^/^<i>_,g^/2<i>_) = (<i)_,g<i)_) = e^. 

Therefore Ej^ < E_ follows. Next we show that Ej^ < E_ by a contradiction. Assume 
that E+ = E_. Then by 1^ we have E+ = {Q^/^^_,Q^/^^_), which implies that 
$_ is a ground state of Q+. In other words, $_ is an eigenvector of Q with eigenvalue 

Ml 



Ea.. Let u 



U2 



be an eigenvector of Q with eigenvalue A: 
A ( ----r + -x^ 1 M + J ( x-^ + -\u = \u. 



(3.2) 
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From the eigenvalue equation (13.21) we can directly see that 



idx' 



:U 



X 



d_ 
dx 



2dx^ 



XA' 



X 

Y 



A-'J 



X 



d 
dx 



■^"4"+{^^"-Y 



1 

2 
A-'J 



u 



d 1 
^^ + 2 



u 



in the sense of distribution. Since x'^Uj.d^Uj/dx'^ G L^(I 

for j = 1,2. By the Sobolev embedding theorem, ui,U2 G C^{ 



A 9 

we see that Uj G W^^^ ( 



Thus $_p,$_p G 



C^(]R) X C^(M) for p = 1, 2. From the fact $_ is odd (resp. $_p is even), $-(0) = 

follows. Therefore $_p satisfies the ordinary differential 



(resp. ((i/(ix) <l>_p(0) 
equations: 





r$-ii 




d 


$-2 




dx 


*'-! 

."^'-2. 








(3.3) 



with the initial condition $„p(0) = and $'„ (0) = 0. Since the right hand side of 
(13.31) is smooth in ($-i, $-2, "^'-i, ^'-2-. t)y the differential equation (13. 3 p has the unique 
solution. Thus $_ = 0, which contradicts ||$-|| = 1. Therefore, E^ < E^. D 

Proof of Theorem \3.1[ Assume that a ^ [5. By Proposition 13. 2[ it is enough to show 
that ker((5 — E) C. l-i^. By Lemma [331 we have £"+ < E_. Hence all the ground 
states are even. Therefore the theorem follows. D 



4. POSITIVITY OF GROUND STATE 



Let 




a„, > 0, 6„ > 0, ?i > 



On > 0,6„ > 0,n > >. 



Then ^^ is a positive cone of H+i and '^q a non-negative cone of H+i. We say that 
^ is non-negative and is denoted by \& > if and only if \E' G "^q ■, and is strictly 
positive and is denoted by \E' > if and only if \& G "io^. A bounded operator T on 
"H+i is positivity preserving if and only if T'^q C ^q , and positivity improving if and 
onlyif T^(+ C ^+. 

Proposition 4.1. Suppose that a hounded self-adjoint operator T is positivity improv- 
ing on l-Lap o-nd ||T|| is an eigenvalue. Then the multiplicity of ||T|| is simple and the 
corresponding eigenvector is strictly positive. 



Proof. See |Far72j . 



D 
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Theorem 4.2. For all t > 0, a = ± and p = 1,2, e~*'^'"p is positivity improving on 
T-Lap. In particular, the lowest eigenvalue of Q^p is simple and corresponding eigen- 
vector is strictly positive. 

Proof. We prove the theorem only for the case of a = + and p = 1. For other cases 
the proof is similar and is left to readers. We shall below show that 6"**^+^ is positivity 
improving. We define 



Ho = A{a*a + ^)T+i, V = ^{aa + a*a*)T+,. 



(4.1) 



Note that Q+i = Hq — V . Since a \n) = ^Jn \n — 1) and a* \n) = y/n + 1 |n + 1), and 
Hq is the multiplication by a{n + |), we see that e~^^° is positivity preserving. Since 



-\4nj 



and 




4n + 2) 


are ana 


e*^ 


'4n)" 



„tv 









c 


An 


+ 2)J 



are analytic vectors oi V, we see that 



P 



^-{aa + 



j=0 

oo 



]An) 




e^+, 



p 



y —{aa + a*a*y 



j=0 





|4n + 2) 



e^+. 



(4.2) 



(4.3) 



From this e*^^o^ '^ ^^ follows. Let ^, $ G '^q ■ By the Trotter-Kato product formula, 
we have 



f^,e-*'3+i$') = hm (^, (e-*^°/^e*^/^ y$) > 0. 

V / i-s>oo ^ 



(4.4) 



Therefore e**^+^ is positivity preserving. Next we show that 6"*^^+^ is positivity improv- 
ing. We can assume that a < P> without loss of generahty. Let P<fc be the projection 
defined by 



P<k 



E,n<k |2^) (4^1 







E4n+2<J4n + 2)(4n + 2| 

It is immediately seen that \I^ > P<fc\& for any \& G '^q^ and e*^/% > (l + tV/j)"^. For 



k' > k, we set v 



\4k) 




and v' 





\4k') 



. Then we have 



v', e 



'v) = lim (v',{e-'^'/^e'^/^yv) > lim (t;', (e-*^«/^P<fc,e*^/^)%) 
> lim fi;',(e-*('='+(^/2))/3/,p^^,gtv/,y \ ^^-t(k'Hi/m ^^ (v',{P<k'{l + tV/j)yv). 

Note that e"*^°(l + tV/j) is still positivity preserving. For all i = 2k' — 2k, we have 
lim {v',{P^k'il + tV/j)yv) > hm {v',jCe{tV/jYv) > lim {v', ,Ce{t{a*f/jYv) 



j-s>oo 



t^ lim .Qj' 



j-s>oo 






V , a 



u 



7i ^™ 



J*^ 



^/ (a*)^k'-ik^ 






v' (a*y''-^''v ] > 0. 
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Thus we have [v', e ^^+^v^ > 0. Similarly 
for all n. Thus e~^^+^ is positivity improving. 



|4n) 






|4n + 2) 



> is derived 
D 



5. No CROSSINGS 

Recall that £'„(«, /3) be the n-th eigenvalue of Q{a,f3), and the c„ : {a,/3) i— )■ 
En{a, /3) G M is an eigenvalue-curve. It is shown that the spectrum of Q is Spec{Q) = 
[Jcr=± p=i 2 '^^^'^(Qpo-)' ^^^ ^11 t^^ eigenvalues in Spec{Qpa) are simple. Now we are 
interested in operators, Q_i — Q+i and Q_2 — Q+2- Let Ao-p(n) be the n-th eigenvector 

1) for n 



I oo 
fn=0 



of Q^p. Then {A^p(n)} 
Theorem 5.1. Assume that 



Spec{Q^p) and A^p(ra) < A^p(?i 



^a/3 > 1 



0,1,2,---. 



1600000000 
Then Q_i — Q+i > A(a, /3) an(i Q_2 — Q+2 > ^(tt, (3), where 

A(a,/3) = 2min{V^, vW"}(v^-l- 1/1600000000) > 0. 
In particular A_i(ra) > A+i(n) + A(a, /3) and A_2(^) > A+2(^) + A(a, f3). 
Proof. We have 



/ 



Q-i ~ Q+i 



\ 



2a 


-70 












-70 




2/3 
-71 


-71 
2a 

-72 


-72 

2/3 

-73 


-73 
2a 

-74 


-74 
2/3 

-75 


-75 
2q 



0\ 



where 7„ = ^^(2^ + 2)(2ra + 3) - ^/{2n + l){2n + 2). We set 

5i = diag[(/3/a)l/^ (a//3)^/^ (/3/a)^/^ («//3)l/^ 
^2 = diag[(«//3)l/^ (/3/a)l/^ («//3)^/^ (/3/a)l/^ 



Then we have 



5l(Q-l - Q+l)Sl = S2iQ-2 - Q+2)S2 

I 



2^/ap -70 
-70 2v^a/3 -71 

—71 2\/ap -72 

—72 2\/a/3 —73 

-73 2v^ 



0\ 







-74 
74 2^0/3 



(5.1) 



(5.2) 



(5.3) 
(5.4) 

(5.5) 
(5.6) 
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We set F = J{{jn)n=o^ 0). Then Si{Q-i-Q+i)Si = 2^0/3- F. Since ^i is self-adjoint 
and invertible, we have 

(g_i - Q+i) > (2v^ - \\F\\)Sf > (2v^ - ||F||) min{ v/^, y^}. 



Similarly we have (Q_2 — Q+2) > {2^/a/3 — \\F\\)m.m{^/aJp, ^J ^jot). Hence it is 
sufficient to prove ||F|| < 2(1 + 1/1600000000). Let v = {vn)^=o e f. Then we have 



I (f , Kv) I = ^{VnlnVfi+l + Vn1nVn+l) 
n=0 

00 00 /■ 2 

< 2^ |t^„|7n|^^n+l| < ^ ian\Vn\'^ + -^ 
n=0 n=0 ^ °" 

for all a„ > 0. So it follows that 



Vn+l\ 



00 2 



n=l 



We split (E?]) as 



A^o 



2 

{v,Kv) I < ao|fo|^ + ^(an + ^^)|WnP 



n=l 



On-l' 



fln-l 



00 2 

E/ I 7n-l\| 1 2 



n=No+l 



fln-l 



with some iVo. We recursively define a„ by 



Oo — 2, a^ 



7n-l 
On-l 



(n = 1,2,3,- ■■ ,iVo), a„ = l (n>A^o + l) 



(5.7) 



(5.8) 



(5.9) 



We can compute the numerical value of a„ from (15. 9|) . e.g. ai = 1.464- ■■, a2 = 
1.305 • • ■ , 03 = 1.228 ■ ■ ■ . We take A'^o = 10000. Then one can easily check that a„ > 
for all n < Nq and a^g > 1. Hence the inequality (15.81) is valid for Nq = 10000 and we 
have 

A'o 00 2 



n=l 

No 



n=No+l 
00 



a-n-1 



<2\vo\' + 2j2\Vn\'+ J2 (l+7n-l)|t'n| 

n=Afo+l 



n=l 



where we used (15. 9p . On the other hand, we have Jn-i = 1 + 
7„_i is monotonously decreasing. Therefore we have 

00 

\{v,Kv)\<{l + jj,^)Y,\Vn\', 

n=0 

which implies that \\K\\ < 1 + 7^^. Note that 



(2n+\/4n2-l) 



^irr-i2- 



In particular 



(5.10) 



7^0 < 7^0-1 < 1 + 
Therefore ||F|| < 2(1 + 1/1600000000). 



(4iVo) 



1 + 



1600000000 



(5.11) 



n 
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The map (a,/3) ^-)■ X(jp{n) = X„p{n,a, (3) is an eigenvalue-curve. It is immediate to 
see the corollary below by Theorem 15.11 . 



;}■ 



Corollary 5.2. Let 

D = {{a,/3) eRxR\a>0,f3>0,ay^ (5, v^ > 1 + 



1600000000- 

Fix p = 1,2. Then two eigenvalue- curves A_p(n) and A+p(n) have no crossing in the 
region D for all n. 

6. Numerical results 

For finite sequences a = (cq, ■ ■ ■ , ctAr-i) and b = {bo, ■ ■ ■ , b]\r), we define the (A^+ 1)- 
dimensional Jacobi matrix, J{a,b), by 

/bo ao \ 

flo bi fli 

ai 62 



J{a, b) 



Vo 



O-N-l b^ 



For cr = ± and p = 1,2, we set a^ = {aij{n))^^Q and 6^ 
finite Jacobi matrix by Q„p{N) = ^J{a^ ,b^p). We set 



^o-p 



(6.1) 



n))^^Q. Define a 



A r/Vwlr../9-lW J"'^^^"^'(2iV + i)'^"'(2iV+^)^ ifiViseven 
+'^ ' 2^ ' ' \xmxi{l3-\2N +l),a-\2N +1)} if iV is odd 



A+2(iV)=A+i(iV) 
1 



A-i(iV) 



(q,/3)^(/3,q) 

Jmin{a-i(2iV + |),/3-i(2Ar+|)} if A^ is even 



2^"" '^ [min{/3-^(2A^ + |),a-i(2A^+|)} if A^ is odd 

A_2(Ar) = A_i(iV) 

(a,/3)^(/3,a) 

and 

|a|a±(A^)| if A^ is even 



^±,i{N) 



\[5\a±{N)\ if A^ is odd. 



^±,2{N) = 5±,i(A^) 



(a,/3)^(/3,a) 



(6.2) 
(6.3) 
(6.4) 
(6.5) 

(6.6) 



Since a/3 > 1, one has Ao-p(A^) = 0{N) — )■ +00 (A^ — )■ +00). Let p„, be the orthogonal 



projection onto e„ = ((5„ 



.iJi=o 



For a self-adjoint operator T, ^niT), n = 1,2, 



denotes the ?7,-th eigenvalue of T counting multiplicity. For n = 0, 1, ■ ■ ■ , A^, we set 

A|,7^^(n) = /i„(g.p(Ar) - (5.p(Ar)p^). 

The eigenvalues of Q^p can be approximated by the eigenvalues of the (A^ -|- 1)- 
dimensional matrix Qcrp{N) in the following sense. 
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Theorem 6.1. Fix N E N, a = ± and p = 1,2. Let n E N be a number such that 

KTn(^) < ^<^p(^)- (6-7) 

Then it follows that 

AS:^(r^) < X.M < KTn(^) (6-8) 

In particular, the error is estimated as \X„p{n) — Xap,N{n))\ < X^^^{n) — X^°pj^{n). 

We give an example below: 

Example 6.2. We set Q± = Q+i{N)±6+i{N)pN- We apply Theorem [6TT] to the case 
a = 1, /3 = 2 and iV = 10. Then A+i(A^) = 5.875 and 

A;;pp^'(0) = 0.366917859 ± 0.000000001, A':^^^^(0) = 0.366917862 ± 0.000000001, 

A:;;pp^'(1) = 2.432911 ± 0.000001, X^+tNi^) = 2.432920 ± 0.000001, 

A^PP^'(2) = 4.7145 ± 0.0001, X^+i%b) = 4.7164 ± 0.0001 

A:;pp^'(3) = 6.2717 ± 0.0001, A'°^'}^(3) = 6.2789 ± 0.0001. 

Since A"pp^(2) < A_|_i(A^) = 5.875, by Theorem 16.11 we have numerical bounds: 
0.36691785 < A^p(O) < 0.36691786, 2.43291 < A^p(l) < 2.43292 and 4.714 < A^p(2) < 
4.717. This example does not include the bound on Ao-p(3), since the condition (16.71) 
is not valid for n = 3. 

Proof of Theorem \6.1[ We prove the theorem only for the case of a = + and p = 1. 
The other cases can be similarly proven. For m, w G £^, we define the operator u Qv : 
P — >■ P by (m © f )$ = (f , $) u, for $ e P. Then operator Q+i can be expressed as 

oo oo 

Q+1 = Q+i(iV) ©0+ ^ 6+i(ra)p„+ ^a+(n)(e„©e„+i + e„+i ©e„). 

n=N+l n=N 

We can show that u Q v + v Q u < eu Q u + e~^v © v for all e > 0. By using this 
inequality, we have 



^ a+(n)(e„ © e„+i + e„+i © e„) < ^ |a+(n)|(e„e„ © e„ + e„^e„+i © e„+i' 

n=N n=N 

oo 

= \a+{N)\eNPN + ^ {en\a+{n)\ + e~Wa+{n - l)\)pn 



n=N+l 



for all en > 0. We take e2n+i = /? and e2n = « for even A^, and e2n+i = « and e2n = P 
for odd A^. Note that |a+(A^)|e7v = 5+i{N). First we consider the case of even A^. 
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Then, we have 

oo 

Y^ {en\a+{n)\ + e~Wa+{n - l)\)pri 

n=N+l 

oo 

= "^if^N+n+ila+iN + n + 1)1 + e^^_^„|a+(A^ + n)\)pN+n+i 

n=0 

oo 

= y^X^N+2n+i\a+{N + 2n + l)\+ e^l^2nk+(^ + 2n)|)pAr+2n.+l 

n=0 

oo 

+ J2i'^N+2n+2\a+{N + 2n + 2)1 + e^l+2n+ll«+(^ + 2?^ + l)\)pN+2n+2 
n=0 

oo 

= ^(/3|a+(iV + 2n + 1)1 + a^'\a+{N + 2n)|)p^+2n+i 

n=0 

oo 

+ ^(a|a+(iV + 2n + 2)1 + /3->+(A^ + 2n + l)\)pN+2n+2- 

n=0 

Since |a+(n)| < 2n + |, we have 

oo 

Y^ (e„|a+(n)| +e-ii|a+(n-l)|)p„ 



n=Af+l 

oo 



< J](^(2iV + 4n + 2 + |) + a-^(2iV + 4n + h)pN+2n+i 



n=0 

oo 



+ J](a(2iV + 4n + 4 + |) + /3-i(2iV + 4n + 2 + |))p;v+2n+2. 



n=0 

By the definition of 6+1(72), we have 
Q+i>Q+iiN)(B0-6+iiN)p^ 



1 °° / 7 3 \ 

+ - ^ ( /3(4iV + 8n + 5) - /3(2iV + 4n + -) - «-i(2iV + 4n + -) j p^+2n+i 

n=0 ^ ^ 

1 °° / 11 7 \ 

+ 2^1 "^^^ + 8n + 9) - a(2Ar + 4n + y ) - r\'2N + 4n + -) j p7v+2n+2 

n=0 ^ ^ 

1 °° 3 

> g+i(iV) © - 6+,{N)pn + - Yi(^ - a~')(2iV + 4n + -)pN+2n+i 

n=0 
1 °° 7 

+ 9 Zl^« - /3"')(2iV + 4n + -)pN+2n+2. 



2 

n=0 



Thus we have Q+i > (Q+i(A^) — 6+i{N)pn) © (A+i(A^)). We can obtain the same 
inequahty for odd A^. In a similar way, we can furthermore obtain the upper bound 
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Q+1 < (,Q+i{N) + 6+i{N)pn) © R{N), where R{N) is an operator such that R{N) > 
A+i(A^). By the min-max principle, we have 

fin{{Q+i{N) - 6+i{N)pn) © A+i(iV)) < /i„(g+i) 
< /i„((g+i(iV) + 6^i{N)p^) © R{N)). 

Suppose that /i„((5+i(iV) + 6+i{N)pn) < A+i(iV). Then 

fin{{Q+i{N) - 6+i{N)p^) © A+i(iV)) = i2n{Q+i{N) - 6+,{N)p^), 
/i„((Q+i(iV) + 5+1 (iVK) © R{N)) = fin{Q+i{N) + 5+1 (iVK). 

This proves (16.81) . D 

7. Concluding remarks 

We can extend non- commutative harmonic oscillators to an infinite dimensional 
version. Let ^ = ©^o-^sym(lR") be the boson Fock space, where Lsyjn(M"), n > 1, 



denotes the set of symmetric square integrable functions, and L^(M'^) = C. Let 
a{f) and a*{f), f G L^(]R), be the annihilation operator and the creation opera- 
tor, respectively, which satisfy canonical commutation relations [a{f),a*{g)] = {f,g), 
[a{f),a{g)] = = [a*{f),a*{g)], and (a(/))* = a*{f). Let dr{u) = J uj{k)a*{k)a{k)dk 
be the second quantization of a real-valued multiplication u. The scalar field is defined 
by (j){f) = 7f(a*(/) + a(/)) and its momentum conjugate by n{f) = ^(a*(/) -a(/)). 
Thus we define the self-adjoint operator 



2' 

on C^ (8) ^. The spectrum of H is not purely discrete. It is interesting to consider the 
existence of a ground state of H and to estimate its multiplicity. 
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